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INTRODUCTION 


Advanced fiber composites are inherently homogeneous material 
whose properties may vary significantly due to small changes in 
fabrication procedure [1]. The nondestructive monitoring of com- 
posite material components for structural integrity assessment is 
important before and during structural service. 

One of the major nondestructive evaluation techniques for mon 
itoring composite materials is ultrasonic testing. It has been 
shown that ultrasonic wave parameters can be correlated with inter 
laminar strength [2], tensile strength [3”5] and transfiber com- 
pression fatigue [1]. Also, knowledge of the wave transmission 
characteristics of composites is important in the interpretation 
of acoustic emission [6] and in dynamic studies such as the assess 
ment of impact damage vulnerability [7]. 

The narrow band ultrasonic longitudinal and shear wave group 
velocities and attenuations in the principal directions of a uni- 
directional Hercules AS/3501-6 graphite fiber epoxy composite and 
in the Hercules 3501-6 epoxy matrix were reported in [8], A 
schematic of the composite laminate is shown in Fig. 1 and the 
examined wave modes are listed in Table 1. The purpose of this 
study is to develop an analysis which relates the composite 
energy loss behavior to the energy loss characteristics of the 
composite constituents. 
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PROPAGATION OF PLANE HARMONIC WAVES IN PRINCIPAL DIRECTIONS OF 
LINEAR VISCOELASTIC SOLID 


If the amplitude of the stress wave Is sufficiently small and 
the wavelength of the harmonic wave is much larger than the dimen- 
sions of the fiber and the fiber spacing, the composite may be 
modeled as a homogeneous anisotropic linear viscoelastic material. 
For plane longitudinal wave propagation along a principal direction 
x in an anisotropic linear viscoelastic rod, there is only one 
nonzero stress component which is defined as a. The constitutive 
equation can be expressed in terms of the complex modulus as [9] 

0 = E(1 + \T])e (1) 


where E and e are the elastic modulus and strain in the x direction, 
respectively, q is the loss factor for longitudinal wave propaga- 
tion in the x direction, and i is the complex number defined by /-T. 
In general, both E and q are frequency-dependent but, in accordance 
with the assumption of linearity, are independent of the amplitude 
of stress and strain. The loss factor is often used to describe 
the energy dissipation property of a material and is defined as 
[ 10 ] 


n = 


aw 

2ttW 


( 2 ) 


where AW is the energy dissipation per unit volume of material per 
cycle and W is the maximum (peak) potential energy within the cycle. 
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For a harmonic longitudinal wave of angular frequency o)(rad/ 

sec) propagating in a linear viscoelastic rod, the phase velocity 

C and the attenuation a are [11] 

P 

c = n k and a = B f- (3) 

P P 

where 

2(1 + and 1 (4) 

/l + n 2 + i ✓! + n 2 + l 

The loss factor n can be expressed explicitly in terms of the 
attenuation and the phase velocity as [12] 

aC 

n = tan (2tan -1 ) . (5) 

Similar procedures can be applied to show that the shear wave 
phase velocity (C p ) s and attenuation a g can be obtained from 
eqns. (3) and (A) by replacing E and n by G and respectively, 
where G is the shear modulus and n g is the corresponding loss 
factor . 

In addition to describing energy dissipation by the loss 
factor r| or the attenuation a, there are several other parameters 
in common use. Relations can be obtained between these various 
parameters and are given in [6]. The loss tangent is simply 

tan 6 = t\ (6) 
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and the logarithmic decrement is 


27raC 

6 ' E. . 

0) 

For small damping, eqns. (6) and (7) can be approximated as 



By assuming an idealized viscoelastic model such as the simple 
Voight model, the damping coefficient C can be expressed as 


C 


In 

03 


The quality factor Q, the critical damping C and the damping 

c 

ratio C are related by 


1 . _ 2C _ ; 

J-2C-5--H 


And, the relaxation time t is defined by 


T a 

E to 


Therefore, the loss factor r) can be related by relatively simpl 
expressions to the other commonly used dissipation parameters. 



LOSS FACTOR OF GRAPHITE FIBER COMPOSITE BASED ON LOSS FACTORS 
OF I TS CONST I TUENTS _______ 

The narrow band group velocity and the attenuation for longi- 
tudinal wave propagation in Hercules 3501-6 epoxy matrix and along 
principal directions of unidirectional Hercules AS/3501-6 graphite 
fiber epoxy composites are reported in [8]. The experimentally 
measured fiber and void volume fractions of the composite were 53% 
and 6 .h% 9 respectively. From those data the loss factor as described 
by eqn. (5) is evaluated for each mode of longitudinal wave propa- 
gation in the composite and the epoxy matrix and are shown in 
Figs. 2 and 3. If the constituent fibers and matrix are assumed 
to be linear isotropic viscoelastic materials, a loss factor of the 
composite for each mode of longitudinal wave propagation can be 
computed, based on the loss factors of the fibers and the matrix. 
Preliminary analysis (similar to that which follows) indicates that 
the resulting relationship between the loss factor of the composite 
and the loss factors of its constituent fibers and matrix is incon- 
sistent with the experimental results. Thus, a new material called 
the “interface material" having unknown properties is introduced as 
an additional constituent of the composite. 

It is postulated that the graphite fiber epoxy composite is a 
combination of fibers, matrix and interface material which, in part, 
can be thought of as regions weakened by voids entrapped during 
fabrication. Simplified models of this combination are shown 
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Figs, 4 and 5. Based on these models, expressions for the composite 
loss factors for different modes of longitudinal wave propagation 
along principal directions will be derived in terms of the loss 
factors of the constituents. 

Loss Factor of Unidirectional Composite for Longitudinal 
Wave Propagation in Fiber Direction 

The equivalent loss factor and elastic modulus of a unidirec- 
tional composite for longitudinal wave propagation in the fiber 
direction are derived next. The composite and its constituent 
fibers, matrix and interface material are assumed to be linear 
viscoelastic with constitutive relations in the fiber direction of 
the composite given by 


E 

E 


E. 


n 


E 

co 


^ ( 1 + i ^ 


(1 + in )e = 

m v m m 


(1 + in. )e. 

in in 


(1 + in )e 

CO CO 


a 

m 


a 

co 


( 12 ) 

(13) 

(14) 

(15) 


where E is the elastic modulus, n is the loss factor, a is the stress 
and the subscripts f, m and in refer to fiber, matrix and interface 
material, respectively, and the subscript co refers to the compos- 
ite in the 0° fiber direction. With regard to Fig. 4, the average 

stress a which acts on the cross-sectional area A of the uni- 
co 

directional composite can be related to the stresses Or , a and a. 

f m in 
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which act on the cross-sectional areas of the fibers A^, the matrix 

A and the interface material A. , respectively. This can be 
m in’ r 7 

written as 


a A = a.Ar + a A + a. A. 
co ft mm in in 


( 16 ) 


By assuming perfect bonding between the fibers, matrix and inter- 
face material, compatibility requires 


e =£,=£. - e 

m f in co 


(17) 


The volume fractions of the fibers, matrix and interface material 
can be written, respectively, as 

A A A 

V„ = 4 - , V = -E- and V. = -ii . (18) 

f A m A in A 

Substituting eqns . (12) to (15) into eqn, (16) gives 


A,. 

E (l + it) ) = E^: ( 1 + iHj:) - !" + E (1 + if) )-r~ + E. (l + it]. ) n 
co co f f A m m A in in A 


i n 


(19) 


Substituting eqn. (18) into eqn. (19) gives 


E (1 + in ) = E,-(l + in JV. + E (1 + in )V + E. (1 + in. )V. 
co co f f f m m m in in in 


( 20 ) 


Separating the real and imaginary parts in eqn. (20) gives 


— LrVr 

CO f f 


m m 


and 


co 


rl r T L. V I 

f f f m m m 


co 


. V. 

in in 

(21) 

v. n. 

in in 

(22) 
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Substituting eqn. (21) into eqn. (22) gives 




E V 
mm 

E f V f 


E. V. 
in in 

E f v f 


n in 


CO 


E V E. V. 

] + m m + in in 


E f v f 


E V 


E f V f 


E V E. V. 

] + m m + in in 


E f v f 


E f v f 


E V E. V. 

I | in m ^ in in 


E f v f 


E f v f 


E. V. 

Assuming « 1 and ' - n . m « 1, eqn. (23) simplifies to 

b f V f Vf 


E V E. V. 

„ x ^ inin„ 

q = + r ' tt ■ q + ■ ■ — q . 

co f E r V £ m E^V^ in 


■f v f 


■f v f 


( 23 ) 


(24) 


Eqns. (21) and either (23) or (24) describe the equivalent elastic 
modulus and the equivalent loss factor of the composite for longi- 
tudinal wave propagation in the fiber direction. If the volume 
fraction of interface material is zero, eqns. (21) and (24) simplify 
to 


E = E-V* + E V 
co f f mm 


E V 

, mm 

q = ru + r ■ q 
co f m 


(25) 

(26) 


Loss Factor of Unidirectional Composite for Longitudinal 
Wave Propagation Perpendicular to Fiber Direction 


For longitudinal wave propagation perpendicular to the fiber 
direction, the equivalent elastic modulus and the equivalent loss 
factor can be obtained from the model shown in Fig. 5« A transverse 
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stress a is applied to the composite. The resulting strains in 

C 9 0 

the fiber, matrix, interface material and composite are, respectively. 


e f E f (t + m f ) 
o 

= C90 

m E m U + "W 


e. 

i n 


C 9 0 

E . (1 + 


i T) . ) 
'in 


in 


a 


C 9 0 




C 9 0 




(27) 

(28) 

(29) 

(30) 


where the subscript C90 refers to the composite in a direction that 
is perpendicular to the fibers and where these constitutive equa- 
tions contain the equilibrium conditions of equal stress on each 
constituent. Because the fibers, matrix and interfaces are in 
series with respect to the loading direction, the composite strain 
is related to the strains in the constituents by 


e = + V e + V. e. 

c 9 o f f mm in in 


(31) 


Substitution of eqns. (27) to (30) into eqn. (31) and separation 
of the real and imaginary parts give 


V. 


1 _ "f . w m + "in (32) 

E C90 U + n l90> " E f<> + E m <> + W E in ( ' + "!7 


and 


n 


v f n f 


v n 

’c 9 o _ f 'f mm ' i n ' i n 

e (i + n* T E,(l + n$) OTTW eTU + n?J 

c 9 o c 9 o t t m m in in 


v. n. 

in in 


. (33) 
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It Is assumed that the damping Is small and so the terms (1 + n?) 
can be approximated by unity; thus eqns. (32) and (33) become 


1 


E 

C 9 0 



(3*0 


T1 

E 


C 9 0 
C 9 0 


V f n f 


v n v. n. 

m m + in in 


m 


E. 


i n 


(35) 


Substitution of eqn. (3*0 into eqn, (35) gives 


n 


E. V 
in m 

E V. 
m i n 


m 


n. 


i n 


C 9 0 


E. V E . V, 
j + m m + inf 


E V. 
m in 


E f V in 


E. V E. V, 
] + i n m + inf 


E. V. 
in f„ 

E,V. n f 
f in 


E V. 
m i n 


E.V. 
f in 


E. V E. V* 
] + m m + inf 


E V. 
m in 


Ej-V . 

f in 


(36) 


Eqns. (3*0 and (36) describe the equivalent elastic modulus and the 
equivalent loss factor of the composite for longitudinal wave prop- 
agation perpendicular to the fiber direction. If the volume fraction 
of interface material V- n is zero, eqns. (3*0 and (36) simplify to 


1 

C 9 0 




E 

v 


E v r 

n - t ) + — m ■ ' n - 
c 9 o m E^V f 
f m 


where the approximation preceding eqn. (2*0 have been used. 


(37) 

(38) 
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RESULTS AND DISCUSSION 


The extensional moduli of the unidirectional AS/3501-6 graph- 
ite fiber epoxy composite as computed from the velocity and density 
measurements are 112 GN/m 2 (16,25 mpsi), 11.67 GN/m 2 (1.69 mpsi) 
and 8.58 GN/m 2 (1.24 mpsi) for the x^ , x 2 and x^ directions, 
respectively. The modulus of the epoxy matrix as computed from the 
velocity and density measurements is 8.94 GN/m 2 (1.29 mpsi). The 
modulus of the fibers is 205 GN/m 2 (29.8 mpsi). The experimental 
frequency-dependent loss factors are given in Figs. 2 and 3- As 
mentioned earlier, the fiber, matrix and void volume fractions are 
53%, 40.6% and 6.4%, respectively; and using eqns. (26) and ( 38 ) 
for a composite with no interface material, the composite loss 
factors could not be verified experimentally. 

The new constituent defined as the interface material is con- 
sidered to be a linear viscoelastic material having, as yet, unknown 
properties. Because the attenuation of the longitudinal wave in 
the x^ direction is much higher than the attenuation in the x 2 
direction, the interface material is assumed to have different prop- 
erties in the and x^ directions. More specifically, the inter- 
face material is assumed to be transversely isotropic with the same 
properties in the x^ and x 2 directions. Therefore, there are seven 
unknown parameters: E* n and q which are transversely isotropic, 

and V. n , V m and . With the experimental data and the additional 
ident i ty 


V, + V 
f m 


+ V. = 1 


( 39 ) 
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eqns. ( 21 ), ( 23 ), ( 3*0 and ( 36 ) can be solved for the values of 
the seven unknowns. Note that in accordance with the transversely 
isotropic assumption for the i nterface mater i a 1 , eqns. ( 3*0 and ( 36 ) 
must be written twice, once each for the and directions. 

Using the results in Figs. 2 and 3 and the seven equations 
cited above, with the nominal composite described by 

E f = 205 §[■ (29.8 mpsi) V f = 53% , V voids = 6 M 

E = 8.9*t ^ (1 .29 mpsi) V = 4 0.6% , (40) 

m m^ m 

the consistent linear viscoelastic model gives 

E f = 205 ^ (29.8 mpsi) V f = 53% 

E = 8.9 1 * 0 .29 mpsi) V = 32.4% 

m m z m 

G N 

3.12 — =- (0.450 mpsi), in x, and x_ directions, 

l z 

1.88 (0.274 mpsi), in x~ direction 

m z 3 

V. = 14. 6% . (41) 

i n 

For consistency, E^., V^. and E m were maintained fixed when going 
from the nominal composite to the consistent linear viscoelastic 
composite. The corresponding loss factors of the interface mate- 
rial in the x^ , x^ and x^ directions are shown in Fig. 6. 

The results given in eqn . (4l ) and Fig. 6 support the concept 
of the interface material as one of the composite constituents. 
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It is important to note that the interface material is an equiva- 
lent volume of material which, if included as a constituent of the 
composite, results in the measured velocities and attenuations. 

Thus, the losses which are attributed to the interface material may, 
in fact, be due to mechanisms such as scattering due to voids and 
delaminations. Although there has been no attempt to relate the 
loss factor of the composite to loss factors of the constituents 
in the shear modes, the analysis developed is still valid for shear 
wave propagation if the extensional variables are replaced by their 
corresponding shear variables. 

For longitudinal wave propagation in the composite, the con- 
tribution from loss factors of the fiber, matrix and interface 
material are weighted as 


A .T) 
m,j m 


where the parameters A r . , A . and A. . are the weighting con- 

f,J m, j i n , j 

stants for the fiber, matrix and interface material, respectively. 
The j subscript distinguishes between weighting parameters for the 
Xj , x^ and x^ directions. According to eqn. (23), the weighting 
parameters for the longitudinal wave in the fiber direction (xj) are 




(^ 3 ) 
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E V 
m m 

E f v f 


m, o 


E V E. V. 

. , mm, in in 
1 + + 


E f v f 


E f v f 


A,. 


in, o 


E in V in 

E f V f 


E V E. V. 

1 + m m . 'n in 

E f V f E f V f 


(44) 


(45) 


For longitudinal wave propagating perpendicular to the fiber direc- 
tion (x 2 or x^), the weighting parameters from eqn. (36) are 


E. Vr 
in f 

E.V. 
f i n 


T , 90 


E. V, E. V 
] + >n f + i n m 


E f V in 


E V. 
m i n 


m, 90 


E. V 
in m 

E,V. 
f i n 


E . Vj. E. V 
j + i n f + i n m 


E,V. 
f i n 


E V. 
m i n 


(46) 


(47) 



90 


] 

E. V, E. V 
in f in m 

+ E.V. E V. 
fin mm 


(43) 


The numerical values for the longitudinal wave weighting parameters 
are summarized in Table 2. 
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The fiber loss factor has been calculated and its behavior 
follows the curve in Fig. 2 within 5% to 10%. So, from Table 2, 
the results indicate that for longitudinal wave propagation in the 
fiber direction, the composite energy dissipation is due mainly to 
the fiber. Also, the composite energy dissipation along the other 
two directions is due mainly to the matrix and the interface 
material . 
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CONCLUSIONS 


Ultrasonic attenuation and velocity measurements were pre- 
viously conducted for longitudinal wave propagation in the Hercules 
3501-6 epoxy matrix and in the principal directions of the uni- 
directional Hercules AS/3501-6 graphite fiber epoxy composite [8]. 

In the present study, the wave propagation results have been inter- 
preted by considering the composite and the epoxy matrix as linear 
viscoelastic media. The loss factors were computed and plotted for 
longitudinal wave propagation in the principal directions of the 
composite and in the epoxy matrix. 

It was observed that the experimental values of the loss 
factors of the composite could not be rationalized to simple paral- 
lel and series models if the composite was assumed to be comprised 
exclusively of fibers, matrix and voids. It was concluded that in 
addition to the fibers and the matrix, a new interface material was 
needed to obtain consistent results between the simple parallel 
and series models and the experiments. The composite, fibers, 
matrix and interface material were each assumed to be linear visco- 
elastic materials. The fibers and the matrix were assumed to be 
isotropic and the interface material was assumed to be transversely 
isotropic. The properties of the interface material were determined 
for longitudinal wave propagation along the principal directions of 
the composite. For longitudinal wave propagation in the composite, 
the contributions of the loss factors of the fibers, matrix and 
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interface material to the composite loss factors were expressed in 
a single equation with weighting coefficients. 

It is important to note that the numerical values of the loss 
factors apply only to the AS/3501 -6 unidirectional composite with 
a 53% fiber volume fraction and a void volume fraction. 

However, the concept of the equivalent dissipative interface mate- 
rial is useful for the quality control of otherwise like composites 
and for consistent wave propagation studies. 
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TABLE I: Propagation and Particle Motion Directions for 

Measurements of Group Velocity and Attenuation [8] 



Propagation 

Particle Motion 



D i rection 

Direction 

Notation 

ro 

n 

x i 

x i 

LW (x j ) 

EH 

X 2 

X 2 

LW(x 2 ) 

WESm 




Lor 

x 3 

x 3 

LW(x 3 ) 



x 2 

SW (Xj ) ; x 2 


X 1 

x 3 

SW(xj); x 3 

<D 



SW(x 2 ); x, 


X 2 

X I 





fO 

<D 

_c 

x 2 

x 3 

SW(x 2 ) ; Xj 

to 

x 3 

X 2 

SW(x 3 ) ; x 2 


x 3 

X 1 

SW(x 3 ); X] 
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TABLE 2: Weighting Factors for Contributions from 

Loss Factors of Composite Constituents 
Towards Loss Factor of Composite 


Direction of 
Long i tud i na 1 
Wave Propagation 

F i ber 

Matrix 

Interface 

Material 

x i 

0.9701 

0.0258 

0.0041 

X 2 

0.0302 

0.4233 

0.5465 

X 3 

0.0222 

0.3111 

0.6667 
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3 Loss factor for longitudinal wave propagation in AS/3501^-6 graphite fiber 
epoxy composite in the indicated directions and 3501-6 epoxy matrix. 






Fig. k Schematic of composite volume element containing 
fiber, matrix and interface material with loading 
in f i ber d i rect ion . 


2k 


"c90 



Fig, 5 Schematic of composite volume element containing 
fiber, matrix and interface material with loading 
perpendicular to fiber direction. 
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Interface Material Loss Factor 


V in (LW(Xj)J 


7] ]u (lw(x, ) a lw(x 2 )) 



2 4 6 8 IxlO 6 2 

Frequency (Hz) 

4 

6 8 IxlO 7 

Fig . 6 

Loss factor of AS/3501-6 Interface material for 
in the indicated directions. 

long i tud i na I 

wave propagation 
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